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ABSTRACT 


During  the  reporting  period  short-period  data  from  several  sites 
have  been  used  to  study  the  nature  of  the  spectra  from  presumed 
explosions,  to  test  a  potential  new  discriminant,  and  to  study 
depth  phases.  Surface-  and  body-wave  magnitude  data  have  been 
obtained  and  used  to  study  regionalization  phenomena.  Prelim¬ 
inary  studies  of  ultra-long-period  data  have  been  undertaken. 
Modifications  of  continental  array  processing  methods  have  been 
completed  and  signal  equalization  studies  initiated.  Upgrading 
of  software  and  hardware  facilities  has  continued. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 
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SUMMARY 


This  is  the  twelfth  Semiannual  Summary  of  Lincoln  Laboratory's  work  for  the  Advanced 
Research  Projects  Agency  on  the  seismic  discrimination  problem  (Vela  Uniform). 

The  Laboratory  has  continued  its  studies  of  short-period  data  and  discriminants.  Emphasis 
has  been  upon  the  simultaneous  utilization  of  data  from  several  scattered  sites.  Cepstral  analysis 
methods  have  been  applied  to  four  Soviet  presumed  underground  explosions  as  recorded  at  Norway, 
LASA,  and  three  United  Kingdom  array  sites.  Reflections  from  the  earth's  surface  above  the 
shot  points  may  have  been  detected.  The  same  data  have  been  used  to  detect  changes  as  a  func¬ 
tion  of  yield  in  short-period  spectra  which  are  consistent  with  those  predicted  by  theory  for  ex¬ 
plosions  (Sec.  1-A).  LRSM  data  from  stations  in  North  America  have  been  used  for  a  preliminary 
test  of  a  discriminant  based  on  body-wave  magnitude,  period,  and  the  minimum  complexity  ob¬ 
served  for  the  set  of  stations  utilized  (Sec.  1-B).  A  limited  study  of  depth  determination  using 
possible  depth  phases  at  LASA  and  NORSAR  has  been  completed.  Utilization  of  two  or  more 
stations  reduces  the  probability  of  erroneously  locating  an  explosion  at  depth  based  on  apparent 
depth  phases  (Sec.  I-C). 

Investigation  of  discrimination  based  upon  surface-  and  body-wave  magnitude  (Mg  —  m^)  has 
continued  with  emphasis  upon  the  regional  effects  due  to  both  receiver  and  source  locations.  A 
study  of  western  United  States  events  using  LRSM  data  verifies  that  presumed  explosions  and 
earthquakes  tend  to  be  separated  by  —  m^  but  that  the  two  populations  are  close  together 
(Sec.  I-D).  Additional  M^  —  m^  data  for  Asia  have  been  measured  at  Norway  and  compared  to 
similar  LASA  data.  A  tentative  interpretation  of  Kurile  Islands  data  in  terms  of  a  downthrusting 
lithosphere  has  been  undertaken  and  preliminary  analysis  indicates  this  will  be  useful  for  region¬ 
alization  studies  (Sec.  I- F). 

Some  studies  of  explosion  characteristics  at  ultra-long  periods  have  been  undertaken. 
Frequency-wavenumber  analyses  of  data  from  NTS  shots  recorded  at  LASA  have  been  computed 
in  the  frequency  range  from  0.01  to  0.05  Hz.  Significant  coherent  signals  were  not  detected  at 
about  0.02  Hz  but  were  detected  at  lower  frequencies  (See.  1-G).  Spectra  have  also  been  obtained 
for  teleseismic  explosions  detected  by  ultra-long-period  instruments  operated  by  the  M.I.T. 

Earth  and  Planetary  Sciences  Department  at  Harvard,  Massachusetts  (Sec.  I-H). 

Array  processing  research  has  included  studies  of  continental- sized  arrays  and  data  equal¬ 
ization  for  array  processing.  Continental  array  beamforming  has  been  modified  to  compensate 
for  variations  in  signal  power  in  time  and  at  different  sites  (See.  Il-A).  Equalization  filters  for 
LASA  have  been  investigated  as  a  preliminary  to  studies  of  deconvolution  filters  which  might 
be  used  to  improve  secondary  phase  detection  (Sec.  Il-B). 

Previously  developed  tools  have  been  applied  to  problems  of  geophysical  interest.  The 
potential  value  of  LASA  and  array  methods  for  studies  of  core  phases  has  been  indicated 
(Sec.  111-A).  The  high-resolution  frequency- wavenumber  method  has  been  applied  to  micro¬ 
barograph  data  (Sec.  Ill- B). 

Sections  III— C  and  111-D  describe  changes  in  our  PDP-7  computer  facility  in  Cambridge  and 
support  software  for  it  as  well  as  for  the  M.I.T.  Information  Processing  System  IBM  360  computers. 

R.  T.  Lacoss 
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GLOSSARY 

ED  R 

Earthquake  Data  Report 

F-K 

Frequency- Wavenumber 

I  PC 

Information  Processing  Center 

LAMA 

Large  Aperture  Microbarograph  Array 

LASA 

Large  Aperture  Seismic  Array 

LRSM 

Long  Range  Seismic  Measurement 

NTS 

Nevada  Test  Site 

OONY 

Oslo  Norway 

SATS 

Semiannual  Technical  Summary 

UKAEA 

United  Kingdom  Atomic  Energy  Authority 

USCGS 

United  States  Coast  and  Geodetic  Survey 

VESPA 

Velocity  Spectral  Analysis 
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SEISMIC  DISCRIMINATION 


I.  SEISMIC  IDENTIFICATION  RESEARCH 

A.  SHORT-PERIOD  SEISMIC  SPECTRUM  OF  EXPLOSIONS 

A  study  has  been  made  of  the  short-period  P-wave  spectra  recorded  at  five  arrays  from 
four  presumed  Soviet  nuclear  tests.  The  purpose  of  the  study  is  to  determine  if  certain  infor¬ 
mation  regarding  the  explosive  source  parameters  can  be  obtained  from  short-period  spectra. 
Specifically,  a  modulation  of  the  spectrum  due  to  surface  reflections  near  the  source  and  a  shift 
in  spectral  content  with  yield  have  been  sought. 

The  data  have  been  taken  from  LASA;  the  United  Kingdom  arrays  at  Yellowknife  Canada, 
Warramunga,  Australia,  and  Gauribidanur,  India;  and  an  array  at  the  LRSM  station  OONY  at 
Lillehammer  near  Oslo,  Norway.  The  global  location  of  these  arrays  and  great  circle  paths  to 
the  source  region  are  shown  in  Fig.  1-1.  The  United  Kingdom  data  were  kindly  supplied  in  digital 
form  by  the  UKAEA  Data  Analysis  ("enter  for  Seismology  at  Blacknest,  England.  It  was  refor¬ 
matted  at  Lincoln  laboratory  in  order  to  make  it  compatible  with  established  analysis  procedures. 

The  LRSM  data  was  digitized  at  Lincoln  Laboratory.  The  data  analysis  console,  described  in 

\ 

a  previous  SATS,  was  used  as  the  primary  analysis  tool. 

It  is  hoped  that  the  use  of  array  beams  in  these  spectral  source  studies  has  minimized  re¬ 
cording  site  crustal  effects  and  enhanced  spectral  features  due  to  the  sources.  Of  course,  this 
assumes  the  crust  will  vary  beneath  the  array  and  be  averaged  out  in  the  beam.  Under  this 
assumption,  the  greater  attenuation  of  recording  site  effects  is  anticipated  as  the  maximum  array 
dimension  increases.  The  maximum  dimensions  of  LASA,  the  United  Kingdom  arrays,  and  the 
Oslo  array  are  roughly  ZOO  km,  25  km,  and  5  km,  respectively.  The  array  beams  were  formed 
by  aligning  the  first  discernible  motion  using  the  data  analysis  console.  A  10- sec  data  time 
window  was  used  for  all  the  spectra  computed.  The  four  presumed  explosions  were  located  by 
LASA  to  lie  in  eastern  Kazakh  within  a  circle  of  2.5°  radius  and  were  assigned  body-wave  mag¬ 
nitudes  of  5.3,  5.4,  5.6,  and  6.1  by  LASA.  The  USCGS  located  the  three  latter  events  to  lie 
within  a  circle  of  0.5°  radius  and  assigned  body  wave  magnitudes  of  5.3,  5.3  and  5.7.  Here,  the 

events  will  be  referred  to  by  their  LASA  magnitudes. 

2 

(‘epstral  analysis  was  performed  on  the  spectrum  of  each  event  as  recorded  at  each  of  the 
array  sites.  The  log  spectrum  was  computed,  a  second  order  trend  removed,  a  taper  to  zero 
at  0.0  and  3.0  Hz  applied,  and  the  spectrum  of  the  log  spectrum  (cepstrum)  computed.  Under  the 
assumption  of  a  plane  wave  elastic  reflection  (pP)  at  the  surface  above  the  source,  the  delay  (r) 
of  this  reflection  is  calculable  from  the  cepstrum.  The  delay  is  related  to  depth  of  burial  and 
other  parameters  by  r  -  2h / (cv  cosi)  where  h  is  the  depth  of  burial,  a  the  comprcssional  wave 
velocity  in  the  overburden,  and  i  the  angle  of  incidence  at  the  receiver.  At  the  epiccntral  dis¬ 
tances  considered  here,  cosi  varies  only  from  about  0.88  to  0.96.  Table  1  lists  the  values  of  r 
to  the  nearest  0.1  sec  which  were  inferred  from  maxima  in  the  computed  cepstra.  In  cases 
where  two  distinct  maxima  were  observed,  two  values  for  r  are  listed. 
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TABLE  1 

DELAYS  INFERRED  FROM  MAXIMA  IN  COMPUTED  CEPSTRA 

Delay  t  (sec) 

LASA  mb 

Norway 

Australia 

India 

Canada 

LASA 

5.3 

No  data 

0.4 

0.6 

0.4 

0.5 

5.4 

0.6 

0.5 

0.4,  0.8 

0.7 

0.4,  0.9 

5.6 

0.6 

0.4 

No  data 

0.7 

0.5 

6.1 

0.5 

0.4,  0.9 

0.5 

0.6 

0.5 

The  cepstral  results  do  not  appear  to  be  definitive.  Although  the  values  obtained  are  rea¬ 
sonable,  the  variation  between  sites  is  greater  than  can  be  accounted  for  by  the  variation  of  the 
angle  of  incidence.  Also,  assuming  the  overburden  velocity  to  be  unchanged,  an  increase  of  r 
with  magnitude  might  be  expected,  since  increased  magnitude  usually  means  increased  depth, 
but  is  not  evident  from  these  data.  Finally,  at  Norway,  Australia,  and  LASA  the  values  of  t 
seem  to  remain  rather  constant,  indicating  a  possible  array  site  effect,  while  at  Canada  and 
India  such  eonstancy  is  not  observed.  This  work  is  continuing  and  will  inelude  analysis  of  earth¬ 
quakes  as  well  as  presumed  explosions. 

Possible  shifts  in  the  spectral  content  of  presumed  explosions  as  a  function  of  yield  have 

been  investigated  by  using  the  ratios  of  short-period  spectra  for  different  events  recorded  at 

3 

each  array.  Displacement  spectra,  based  on  Haskell's  analytical  approximation  to  experimen¬ 
tal  curves  for  reduced  displacement  potentials  of  contained  underground  explosions,  were  used 
as  a  basis  for  comparison.  Haskell’s  displacement  spectra  for  granite  at  various  yields  are 
shown  in  Fig.  1-2.  Note  that  there  is  considerable  change  in  the  shape  of  the  short-period  spec¬ 
trum  as  a  function  of  yield. 

In  Fig.  1-3  the  log  spectral  ratios  log  A(f)^  ^/A(f)j.  4  and  log  A(f)^  ^/A(f)^  ^  computed  for 
LASA  and  Australia  arc  plotted  over  the  bands  0.6  to  2.0  Hz.  Here  A(f)  is  the  speetrum  of  the 
seismometer  output  whieh  is  roughly  proportional  to  ground  velocity  in  this  frequency  band.  The 
subscripts  indicate  the  LASA  magnitude  of  the  event.  By  taking  spectral  ratios  at  a  fixed  re¬ 
cording  site  of  events  from  similar  locations,  cancellation  of  the  spectral  effects  of  the  seis¬ 
mometer  and  the  earth  transmission  should  take  plaee,  leaving  only  those  due  to  variation  of  the 
source  spectrum.  The  solid  straight  line  is  a  least  squares  fit  to  the  observations.  The  dashed 
line  represents  the  spectral  ratio  based  on  Haskell's  expressions  for  granite  and  the  magnitude 
yield  relation  =  3.7  f  log  (kilotons).  The  change  in  slope  of  the  fitted  line  with  increasing 
magnitude  differential  is  apparent.  Tabic  11  summarizes  the  spectral  ratio  data  for  all  of  the 
array  sites.  The  slope  of  the  line  fitted  to  the  observed  spectral  ratios  is  given  for  each  array 
and  values  derived  from  Haskell's  expressions  are  given  for  comparison.  At  each  array  (where 
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SLOPE 

TABLE  II 

OF  LEAST  SQUARES  FIT  TO  RATIO  OF 

SPECTRA 

LASA 

Ratio 

Norway 

Australia 

India 

Canada 

LASA 

Haskel  1 
Theoretical 
Slope 

5. 6/5. 4 

0.006 

-0.201 

Na  data 

+0.011 

+0.023 

-0.083 

6. 1/5.4 

-0.229 

-0.398 

-0.280 

-0.119 

-0.310 

-0.367 

data  is  available)  the  slope  of  the  fitted  line  decreases  with  increasing  magnitude  differential. 

This  is  taken  as  a  good  indication  of  the  existence  of  source-related  information  in  the  short- 
period  spectrum  at  telescismic  distances. 

J.  Fils  on 

B.  DISCRIMINATION  USING  WORLDWIDE  OBSERVATIONS  OF  PERIOD, 

COMPLEXITY  AND  BODY-WAVE  MAGNITUDE 

It  has  been  found,  using  LASA  data,  that  valuable  discriminants  can  be  based  on  the  use  of 
the  spectral  ratio  and  the  relation  between  surface- wave  and  body- wave  magnitude,  Mg  —  m^. 
However,  these  discriminants  are  not  useful  at  LASA  below  a  certain  magnitude  threshold.  It 
is  not  clear  that  the  spectral  ratio  defined  for  LASA  data  will  provide  effective  discrimination 
results  when  obtained  at  other  stations.  This  statement  is  based  on  preliminary  data  from 
NOKSAK.  It  is  possible  to  use  the  Mg  —  discriminant  effectively  at  other  stations.  However, 
the  utility  of  this  discriminant  is  limited  by  the  difficulty  in  detecting  the  Rayleigh  surface  wave 
due  to  the  poor  signal-to-noise  ratio  on  the  long-period  instruments,  as  compared  to  that  on  the 
short-period  sensors.  Thus,  in  attempting  to  lower  the  identification  threshold  by  using  world¬ 
wide  observations,  it  is  desirable  to  utilize  discriminants  other  than  those  mentioned  previously. 

One  effective  way  of  lowering  the  identification  level  is  to  use  discriminants  based  solely 
on  the  gross  characteristics  of  the  P-wave,  due  to  the  high  signal-to-noise  ratio  on  t lie  short- 
period  seismometers.  A  preliminary  investigation  of  such  a  discriminant  will  now  be  described. 
The  quantity  P  is  defined  as 

p  =  cmin(T-°-5)2 

where  C ^ ^  is  the  minimum  complexity  observed  across  a  network,  the  complexity  at  the 
station  is  defined  as 

/235l»jMdt 

fo  |Dj(t)|di 

D  (t)  is  the  P-wave  data  at  the  ^  station,  t  is  the  time  in  seconds,  the  onset  time  of  the  P-wave 
corresponds  to  zero  time,  T  is  the  average  of  the  dominant  periods  observed  across  the  network 
and  the  dominant  period  is  the  period  associated  with  the  largest  peak- to- peak  excursion  of  the 
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P-wave  amplitude  in  the  first  few  cycles  following  the  onset  time.  The  actual  discrimination  is 
performed  by  plotting  P  vs  m^,  where  fh^  is  the  average  of  the  body- wave  magnitudes  observed 
across  the  network  and  is  defined  in  the  usual  way  for  teleseismic  events. 

The  results  of  applying  this  discriminant  to  ten  earthquakes  and  eight  presumed  underground 
nuclear  explosions,  all  from  the  Central  Asian  region,  using  film  data  recorded  at  LRSM  sta¬ 
tions  in  North  America,  are  shown  in  Fig.  1-4.  An  average  of  about  six  stations  was  used  for 
each  observation  point.  These  preliminary  results  are  encouraging  since  there  is  no  overlap 
between  the  two  populations.  However,  it  is  desirable  to  perform  the  experiment  using  more 
events,  and,  in  particular,  events  with  body- wave  magnitudes  in  the  3.5  to  4.5  range.  It  is  im¬ 
portant  to  note  that  the  discriminant  based  on  P  vs  m^  may  be  measured  visually  by  an  operator 
using  film  data  and  does  not  necessarily  require  elaborate  digital  or  analog  signal  processing 
equipment.  This  may  be  an  advantage  if  only  film  records  from  a  network  of  stations  are  avail¬ 
able.  However,  computer  compatible  recordings  would  significantly  reduce  the  labor  required 
for  complexity  measurements  and  would  probably  increase  accuracy  and  consistency. 

J.  Capon 
L.  Lande 


C.  LASA-NORSAR  DEPTH  PHASE  STUDY 

If  it  can  be  shown  that  an  event  has  occurred  at  a  depth  greater  than  a  few  kilometers,  such 
an  event  can  be  identified  as  an  earthquake.  The  major  difficulty  with  this  method  is  the  possi¬ 
bility  of  erroneously  locating  an  explosion  at  a  depth  greater  than  a  few  kilometers.  This  prob- 

4 

lem  was  graphically  demonstrated  in  a  recent  report  where  a  presumed  explosion  clearly  showed 
an  arrival  at  LASA  which  could  have  been  interpreted  as  pP  or  sP  and  would  have  located  the 
event  deep  in  the  crust.  The  usual  suggestion  is  that  such  problems  can  be  circumvented  by 
requiring  corroborating  apparent  depth  phases  at  several  other  receiver  sites.  We  have  now 
completed  a  limited  study  to  evaluate  this  suggestion  using  only  LASA  and  NOKSAR  data. 

A  total  of  41  earthquakes  ana  six  presumed  explosions  have  been  found  with  usable  short- 

4 

period  NORSAR  data,  LASA  beams  had  previously  been  formed  and  possible  depth  phases  iden¬ 
tified.  These  events  were  all  from  eastern  Europe  and  Asia.  Beams  were  formed  using  the 
eleven  available  sensors  at  NORSAR  and  an  analyst  identified  possible  depth  phases  on  those 
beams.  The  presumed  explosion  which  showed  a  possible  depth  phase  at  LASA  was  not  in  this 
population.  However,  a  different  presumed  explosion  showed  a  possible  depth  phase  at  NORSAR 
which  was  not  corroborated  at  LASA. 

Table  111  is  a  summary  of  depths  implied  by  depth  phase  picks  for  the  37  earthquakes  for 
which  either  site  showed  a  provisional  depth  phase.  The  depth  shown  is  that  which  follows  from 
assuming  that  the  phase  was  pP.  If  the  phase  were  sP,  the  depth  would  be  reduced  by  a  factor 
of  roughly  2/3.  A  zero  entry  indicates  no  depth  phase  was  picked  by  the  analyst.  Events  in  the 
table  have  been  identified  as  showing  good,  fair,  or  poor  agreement  between  the  LASA  and 
NORSAR  depths.  Good  or  fair  agreement  was  obtained  for  54  percent  of  the  events.  Compari¬ 
son  with  USCGS  events  indicate  that  some  events  with  poor  agreement  probably  had  a  correct 
depth  phase  picked  at  one  or  both  sites. 
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DEPTH  DATA  FOR  USCGS 

TABLE  III 

AND  LASA-NORSAR  BASED  UPON  POSSIBLE  pP  PICKS 

Events  with  Good 

Events  with  Foir 

Events  with  Poor 

LASA- 

NORSAR  Agreement 

LASA- 

NORSAR  Agreement 

LASA- 

NORSAR  Agreement 

LASA 

NORSAR 

USCGS 

LASA 

NORSAR 

USCGS 

LASA 

NORSAR 

USCGS 

150 

160 

155 

34 

40 

50 

0 

90 

- 

580 

580 

582 

30 

55 

33R 

15 

90 

- 

450 

470 

463 

23 

30 

140 

40 

161 

510 

510 

511 

15 

20 

0 

8 

- 

20 

25 

21 

55 

30 

0 

15 

- 

120 

110 

118 

30 

50 

150 

15 

134 

15 

15 

33R 

50 

33 

0 

15 

20 

170 

170 

160 

17 

11 

110 

0 

113 

30 

25 

- 

0 

80 

16 

20 

17 

- 

20 

100 

33R 

46 

44 

- 

33 

10 

50 

45 

50 

- 

115 

15 

10 

50 

0 

- 

75 

0 

94 

0 

21 

33R 

70 

20 

33  R 

0 

160 

- 

Extrapolating  to  a  system  using  many  sites,  it  appears  that  50  percent  or  more  of  the  earth¬ 
quakes  might  be  identified  by  depth  phases  without  erroneously  locating  any  explosions  at  depth 
This  could  be  very  significant  since  many  earthquakes  identified  in  this  way  might  be  trouble¬ 
some  events  for  other  criteria.  An  experiment  using  many  sites  and  many  events  is  required. 

R.  T.  Lacoss 
L.  Lande 

D.  DISCRIMINATION  OF  NORTH  AMERICAN  EVENTS  USING  BODY- 
AND  SURFACE-WAVE  MAGNITUDES 

A  valuable  discriminant  for  distinguishing  between  earthquakes  and  underground  nuclear 

explosions  is  based  on  the  relation  between  the  surface-wave  mangitude  M  .  and  the  body-wave 

4  5 

magnitude  m^.  It  has  been  found  1  that,  above  a  body-wave  magnitude  of  about  4.9  and  for  pre¬ 
sumed  Soviet  explosions,  there  is  near-perfect  separation  between  the  two  populations  using  the 


5 


Section  I 


discriminant  oven  if  data  from  only  a  single  site  (LASA)  is  utilized.  However,  recently  it  has 
been  proposed  by  Liebermann  and  Pomeroy^  that  this  discriminant  will  not  provide  a  perfect 
separation  between  the  two  populations,  particularly  in  the  body-wave  magnitude  range  m^  <  4. 

In  order  to  investigate  this  problem,  measurements  of  Mg  and  m^  were  made  on  24  earth¬ 
quakes  and  nine  presumed  underground  nuclear  explosions  in  North  America,  using  film  data 
recorded  at  LRSM  stations  in  North  America.  The  average  value  Mg  of  the  surface-wave  mag¬ 
nitudes  observed  across  the  network  was  used,  as  well  as  the  average  body- wave  magnitude  m^ 
The  value  of  m^  at  each  station  was  computed  according  to  the  procedure  suggested  by  Evernden7 
for  near-regional  distances  in  the  United  States.  An  average  of  about  5  LRSM  stations,  varying 
in  distance  from  about  2°  to  40°  from  the  epicenter  of  the  event,  were  used  for  each  measure¬ 
ment  of  and  m^.  The  results  of  the  measurement  are  shown  in  Fig.  1-5,  along  with  the 
Gutenberg-Riehter  relationship,  Mg  =  t.59m^  —  3.97.  It  is  observed  from  Fig.  1-5  that  there  is 
not  a  perfect  separation  between  the  two  populations.  The  open  circle  which  lies  above  the 
Gutenberg-Riehter  line  corresponds  to  the  Gasbuggy  explosion.  The  two  populations  appear  to 
be  very  close  together,  particularly  at  the  lower  body-wave  magnitudes.  They  are  much  closer 

5 

together  than  similar  populations  for  central  Asia.  This  result  is  in  reasonably  good  agreement 
with  that  given  by  Liebermann  and  Pomeroy/5  It  is  possible  that  the  two  populations  may  be  close 
together  due  to  the  triggering  of  earthquakes  in  the  vicinity  of  the  Nevada  Test  Site  by  the  under¬ 
ground  nuclear  explosions.  Unfortunately,  our  data  does  not  include  any  very  low  magnitude 

explosions.  T  ^ 

r  J.  Capon 

L.  Lande 

E.  DISCRIMINATION  USING  BODY-  AND  SURFACE-WAVE  MAGNITUDES 
MEASURED  AT  NORSAR 

Mg  and  measurements  have  been  made  at  NORSAR  for  33  shallow  (depth  <  50  km)  earth¬ 
quakes  and  four  presumed  underground  nuclear  explosions  occurring  outside  the  western  United 
States  and  Mexico.  These  are  shown  in  Fig.  1-6.  The  M  value  is  the  average  of  the  Mg  value 
determined  at  each  of  the  three  temporary  long- period  sites  established  in  Norway  by  Lincoln 
Laboratory.  Previously-published  data  used  measured  at  only  one  site,  consisted  of  many 
fewer  events,  and  did  not  include  NORSAR  values. 

Unfortunately,  the  values  could  only  be  determined  exactly  for  two  of  the  four  pre¬ 

sumed  explosions  examined.  The  individual  outputs  of  the  NORSAR  short-period  vertical  instru¬ 
ments  were  clipped  for  the  presumed  explosion  which  is  plotted  at  =  6.0.  We  would  estimate 
the  true  value  of  for  this  event  as  6.2.  The  surface  waves  from  one  of  the  three  other  pre¬ 
sumed  explosions  could  not  be  positively  identified  above  the  noise.  Only  an  upper  bound  on  Mg 
could  be  determined  for  this  event  and  the  event  is  shown  with  an  arrow  pointing  down  to  indi¬ 
cate  this  fact. 

We  now  have  sufficient  data  from  the  Kurile  Islands  and  central  Asia  to  undertake  a  prelim- 

5 

inary  comparison  of  the  NORSAR  data  with  similar  data  obtained  at  LASA.  The  NORSAR  data 
are  presented  in  Figs.  1-7 (a)  and  1-7 (b).  The  solid  straight  line  is  the  Gutenberg-Riehter  rela¬ 
tionship  between  Mg  and  m^.  The  NORSAR  Kurile  Islands  data  cluster  about  mb  =  6  and  the 
central  Asia  earthquake  data  are  clustered  about  =  5.  A  comparison  with  LASA  data  suggests 
some  hypotheses  which  might  be  substantiated  using  considerably  more  data  with  a  wider  variation 
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in  recorded  at  both  sites.  First,  for  the  Kurile  Islands  events  in  the  vicinity  of  m^  -  6,  the 
LASA  M  values  lie  below  the  Gutenberg-Kichter  line,  while  the  NORSAR  vaJ"»s  lie  nearly  on  it. 
The  situation  is  changed  in  central  Asia.  For  this  region,  the  NORSAR  values  still  tend  to  lie 
on  the  Gutenberg-Richter  line  but  the  LASA  values  tend  to  be  above  the  line.  This  shift  at 
NORSAR  brings  the  earthquake  population  and  explosions  population  closer  together. 

The  reasons  for  the  apparent  average  differences  between  NORSAR  and  LASA  M  data 

are  not  entirely  clear  at  this  time.  However,  one  factor  has  been  investigated  using  USCGS  and 
NORSAR  data  for  the  11  Kurile  Islands  earthquakes.  The  mean  of  (NORSAR)  —  (USCGS) 
for  the  11  events  is  0.44.  The  structure  under  the  Kurile  Islands  has  been  portrayed  as  a  low 
attenuation  lithospheric  plate  thrusting  downward  into  a  relatively  high  attenuation  upper  mantle. 
The  seismic  activity  occurs  along  the  upper  edge  and  in  any  graben  features.  A  plausible  expla¬ 
nation  for  the  higher  values  at  NORSAR  is  that  the  body- wave  propagates  through  the  upper 
few  hundred  kilometers  of  the  mantle  in  the  low  attenuation  lithosphere  rather  than  the  surround¬ 
ing  high  attenuation  upper  mantle. 

Figure  I-7(c)  shows  the  NORSAR  Mg  vs  the  USCGS  m^  for  the  Kurile  events.  Comparison 

with  Fig.  1-7 (a)  makes  the  shift  in  m.  values  obvious.  The  m,  values  do  not  cluster  at  m,  =  6.0 

b  b  b 

but  range  over  the  interval  5.0  to  6.0.  The  scatter  also  seems  to  be  significantly  reduced.  Much 
of  the  scatter  in  Fig.  1-7 (a)  may  be  due  to  large  variations  in  body- wave  attenuation  resulting  from 
small  changes  in  epicenter.  This  is  possible  since  small  epicenter  shifts  can  significantly  affect 
the  amount  of  P-wave  travel  in  the  low  attenuation  lithosphere.  The  USCGS  is  not  as  sensi¬ 
tive  to  this  since  it  is  an  average  over  several  stations  at  different  azimuths.  All  of  this  seems 
to  indicate  that  for  a  small  complex  region,  sucli  as  the  Kuriles,  the  difficulty  is  using  M  — 
from  a  single  station  as  a  discriminant  may  result  more  from  scatter  in  m^  values  than  in  M 
values. 

To  ascertain  the  extent  to  which  the  scatter  of  the  NORSAR  central  Asia  events  is  due  to  a 
single  station  determination  of  m^,  we  used  the  multistation  USCGS  m^  for  these  events.  Fig¬ 
ure  I-7(d)  shows  the  result.  Scatter  is  reduced  and  the  separation  between  earthquakes  and  pre¬ 
sumed  explosions  somewhat  enhanced.  However,  the  m^  values  for  USCGS  still  cluster  about 
5.0  as  the  NORSAR  values  did.  Since  the  central  Asia  region  encompasses  a  much  larger 

area  than  the  Kurile  Islands,  the  M  -  m,  relation  for  earthquakes  is  not  as  well  defined  and 

s  b  4  4 

scatter  is  not  reduced  as  significantly.  These  results  agree  with  those  previously  obtained  for 
a  similar  experiment  using  LASA  data.  R  w  Ward  (M  j  T  Earth 

and  Planetary  Sciences  Department) 

F.  SURFACE-  AND  BODY-WAVE  MAGNITUDE  OBSERVATIONS 
REPORTED  TO  USCGS 

A  study  is  under  way  to  determine  the  regional  effects  on  magnitude  obtained  from  surface 
waves  and  body  waves  for  earthquakes.  To  avoid  the  problems  associated  with  obtaining  a  large 
data  base,  such  as  acquiring  seismograms  and  reading  them,  the  data  on  surface-wave  and  body- 
wave  magnitudes  published  by  the  USCGS  in  their  Earthquake  Data  Report  publications  are  being 
used.  For  about  a  year,  the  USCGS  has  been  publishing  surface-wave  magnitudes  associated  with 
approximately  five  percent  of  their  published  epicenters.  Over  the  period  of  one  year,  it  was 
possible  to  obtain  surface-wave  magnitudes  on  over  400  events  and  well  over  a  thousand  individual 
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station  surface-wave  magnitudes.  These  data  have  been  reduced  to  punched  cards  and  presently 
represent  our  initial  data  base.  Periodic  updating  will  occur  and  enlarge  the  data  base  to  an 
expected  several  thousand  readings. 

A  preliminary  analysis  of  data  for  epicenters  in  the  Japan  region  has  produced  some  results. 
When  the  individual  surfaee-wave  magnitudes  are  plotted  versus  the  USCGS  body-wave  mag¬ 
nitudes  m^,  it  can  be  seen  that  stations  in  Europe  have  higher  values  for  the  same  event 
than  stations  in  North  America.  This  effect  is  in  substantial  agreement  with  the  results  obtained 
by  R  Ward  from  Kurile  Islands  earthquakes  recorded  at  NORSAR  and  LASA  and  discussed  in 
Sec.  I-E  of  this  report.  In  addition  to  our  being  able  to  obtain  a  knowledge  of  source  and  receiver 
regional  variations,  we  expect  that  it  will  be  possible  to  determine  station  magnitude  corrections 


for  M  and  m,  . 
s  b 


R  M.  Sheppard 


G.  LONG-  AND  ULTRA- LONG- PERIOD  CHARACTERISTICS  OF  NTS  EVENTS 
RECORDED  AT  LASA 


The  distance  from  NTS  to  Montana  LASA  is  in  the  range  12°  to  15°.  This  is  considerably 
less  than  teleseismic  distance  and  so  very  little  interest  has  previously  been  shown  in  NTS  events 
recorded  at  LASA.  However,  the  application  of  array  processing  techniques  to  surface  waves 
generated  by  nuclear  explosions  at  this  relatively  short  distance  may  indicate  discriminants 
which  could  be  applied  on  a  global  scale.  Furthermore,  the  detection  of  waves  with  periods  up 
to  100  secs  may  be  possible  by  bcamforming. 

9 

Frequency- wavenumber  (F-K)  spectra  have  been  calculated  for  the  vertical  components  of 
the  NTS  events  ZAZA,  LANPHER  and  SLED,  recorded  at  Montana  LASA  The  advantage  of  F-K 
spectra  is  that  the  full  capabilities  of  the  array  to  detect  coherent  propagating  waves  are  utilized. 
Figure  1-8  shows  the  high  resolution  F-K  spectrum  for  ZAZA  at  a  frequency  of  0.05  Hz.  This 
corresponds  to  the  peak  energy  in  the  power  speitrum  of  the  event.  High  resolution  F-K  spectra 
were  calculated  for  several  frequencies  down  to  0.01  Hz  and  it  appears  that  measurable  coherent 
energy  persists  out  to  very  low  frequencies.  Figure  1-9  shows  the  F-K  spectrum  at  a  frequency 
of  0.1  5  Hz.  Tentatively,  it  appears  that  coherent  energy  from  the  events  was  detectable  even  at 
0.01  IIz.  At  such  a  low  frequency  the  resolution  of  LASA  is  poor  and  so  no  great  significance 
can  be  attached  to  the  precise  values  of  azimuth  and  phase  velocity  obtained  at  low  frequencies 
when  compared  with  the  same  quantities  measured  at  higher  frequencies.  An  apparent  loss  of 
coherent  signal  is  encountered  in  the  region  of  0.02  Hz.  Figure  I  - 1  0  show's  the  high  resolution 
F-K  spectrum  at  0.02  Hz  and  there  is  certainly  no  discernible  coherent  energy. 

A  plot  of  average  power  versus  frequency  is  shown  in  Fig.  I- 11  for  ZAZA.  The  average 
power  is  the  power  spectrum  averaged  over  the  individual  sensors  of  LASA  and  has  contributions 
from  both  coherent  and  incoherent  disturbances.  The  values  are  normalized  to  the  power  at  a 
frequency  of  0.05  Hz.  A  dip  in  average  power  is  evident  at  frequencies  around  0.02  IIz.  The 
lower  curve  in  the  same  figure  is  a  measure  of  the  ratio  of  the  coherent  power  to  total  power 
as  determined  by  F-K  analysis  and  a  distinct  drop  in  this  value  is  seen  at  frequencies  around 
0.021  IIz,  corresponding  to  the  notch  in  the  average  power  spectrum.  Figure  1-12  shows  the 
same  measurements  for  the  event  SLED  and  a  very  similar  pattern  emerges. 
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The  above  results  suggest  that  the  LASA  long-period  system  may  be  used  to  study  waves 
with  frequencies  down  to  0.01  Hz,  using  frequency  domain  beamforming  and  F-K  spectra. 

H.  Mack  (M.  1.  T.  Earth 

and  Planetary  Sciences  Department) 

H.  LONG-  AND  ULTRA- LONG- PERIOD  SIGNALS  RECORDED  FROM  EXPLOSIONS 
AT  TELESEISMIC  DISTANCES 

The  M.  I.T.  mercury  tube  tiltmeter  has  been  operating  for  about  one  year  at  the  Agassiz 
Station,  Harvard,  Massachusetts.  A  companion  Geotech  vertical  Press-Ewing  has  been  oper¬ 
ating  for  about  nine  months.  During  this  time,  four  U.  S.  underground  tests  have  been  recorded 
on  the  tiltmeter,  including  one  from  Amchitka  Island.  The  vertical  component  has  recorded  two 
additional  U.  S.  tests,  including  one  from  Colorado,  and  a  Russian  test  from  Novaya  Zemlya.  All 
three  components  have  been  bandpass  filtered  (10-  to  80-second  periods)  to  enhance  long-period 
surface  waves,  so  body  waves  from  small  events  are  not  recorded.  Detection  thresholds,  there¬ 
fore,  are  based  on  recording  surface  waves.  At  35 °,  we  estimate  that  the  tiltmeter  will  record 
Rayleigh  waves  from  underground  explosions  greater  than  500  kt,  whereas  the  vertical  compo¬ 
nent  is  sensitive  enough  to  record  events  down  to  40  kt.  For  earthquakes  at  about  the  same 
distance,  the  vertical  component  has  recorded  surface  waves  from  a  Santa  Rosa,  California 
aftershock  with  M^  =  3.4,  more  than  one  magnitude  unit  less  than  the  minimum  recorded  on  the 
tiltmeter. 

In  their  study  of  long-period  discrimination,  Molnar,  et  al.  defined  amplitude  ratios  for 

Rayleigh  and  Love  as  A^  =  A^  >2/^40-60  anc*  ~  ^17-25^40-60  rcsPectively>  where  A1  ? 

is  the  maximum  amplitude  of  the  Rayleigh  wave  at  period  1  9  to  22  sec,  and  so  on.  They  found 
that  A^  =  6  completely  separated  explosions  from  earthquakes,  but  that  no  separation  could  be 
found  for  A^  .  These  results,  and  the  lack  of  sensitivity  of  the  tiltmeter  compared  to  the  vertical 
long-period  seismometer,  suggest  that  the  tiltmeter  will  be  of  minimal  use  for  discrimination 
at  present.  The  long-period  vertical,  however,  seems  capable  of  providing  a  single  station 
teleseismic  discriminant,  especially  if  installed  at  a  quiet  site.  For  example,  the  instrument 
used  by  Molnar,  et  al.  was  installed  deep  in  a  mine  with  careful  temperature  and  pressure  control 
at  Ogdensburg,  New  Jersey,  and  has  recorded  events  at  30°  to  40°  of  M?  as  low  as  2.7,  better 
than  half  an  order  of  magnitude  below  the  smallest  recorded  at  Harvard.  In  general,  more  test¬ 
ing  of  discriminants  utilizing  ultra-long-period  waves  remains  to  be  done. 

J.  Derr  (M.  I.  T.  Earth 

and  Planetary  Sciences  Department) 
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Fig.  1-1.  Glabal  location  af  arrays  from  which  short-period  data  was  obtained. 


Fig.  1-2.  Displacement  spectra  based  an  Haskell’s 
analytic  approximation  for  various  yields  in  granite. 
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Fig.  1-3.  Lag  ratios  af  presumed  explosion  spectra  ot  (a)  Australio 
and  (b)  LASA.  Solid  straight  lines  are  fitted  ta  data;  dashed  lines 
are  computed  from  Hoskell's  opproximotion . 
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Fig.  1-4.  Discrimination  based  on  worldwide  observa-  Fig.  1-5.  Surface-wave  vs  body-wave  magnitude  far 
tions  of  period,  complexity,  ond  body-wave  mognitude.  North  American  events  using  LRSM  data. 


Fig.  1-6.  Surface-wave  vs  bady-wave  magnitude 
using  NORSAR  data.  Al  I  eventsore  outside  western 
United  States  and  Mexico. 
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Fig.  1-7.  NORSAR  surface-wave  magnitude  vs  bady-wave  magnitude. 

a.  Kurile  Islands  events  with  NORSAR  bady-wave  magnitude. 

b.  Central  Asia  events  with  NORSAR  bady-wave  magnitude. 

c.  Kurile  Islands  events  with  USCGS  bady-wave  magnitude. 

d.  Central  Asia  events  with  USCGS  bady-wave  magnitude. 
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Fig.  1-8.  High  resolution  wavenumber  spectrum  of  vertical  component  at  LASA 
for  NTS  event  ZAZA  at  0.05  Hz.  Contours  ore  in  dB. 
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Fig.  1-9.  High  resolution  wavenumber  spectrum  of  vertical  component  at  LASA 
far  NTS  event  ZAZA  at  0.015  Hz.  Contours  are  in  dB. 


15 


, 


Section  I 


hi-g-aotTl 


^ _ J  J 

\  < 

A  Vc 

:)Jo  c 

lV0>J 

Ni  s 

?  _ _ _ ✓  ' _ ^ 

•  0.015  -0.000  0.015 


WAVENUMBER  (cycles/km) 


Fig.  1-10.  High  resolution  wavenumber  spectro  of  vertical  component  ot  LASA 
for  NTS  event  ZAZA  ot  0.02  Hz.  Contours  ore  in  dB. 
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Fig.  1-1  1  .  Averoge  power  spectrum  and 
coherency  measure  for  NTS  event  ZAZA 
recorded  by  the  LASA  LP  system. 


Fig.  1-12.  Average  power  spectrum  ond 
coherency  meosure  for  NTS  event  SLED 
recorded  by  the  LASA  LP  system. 
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RATIO  OF  COHERENT  POWER 
TO  TOTAL  POWER 


II.  SIGNAL  PROCESSING  RESEARCH 


A.  CONTINENTAL  APERTURE  SEISMIC  ARRAYS 


A  description  of  the  continental  aperture  seismic  array  computer  program  has  been  given 
1 1 

previously.  It  will  be  recalled  that  the  purpose  of  this  program  is  to  determine  the  depth  of 

an  event  by  recognizing  the  arrival  o:  the  pP  phase,  making  use  of  P-pP  differences  in  velocity, 

and  to  determine  the  P-wave  source  structure  of  the  event  by  essentially  steering  many  beams 

in  the  vicinity  of  the  event.  The  data  consist  of  the  digitized  and  merged  records  from  LASA 

and  LRSM  sites  located  in  North  America.  The  results  obtained  for  processing  the  27  October 

1 1 

1966  Novaya  Zemlya  event  with  this  program  were  presented  previously. 

Several  experiments  have  been  performed  recently  to  determine  the  validity  of  the  results 
obtained  using  the  computer  program.  These  results  will  now  be  described.  One  of  the  exper¬ 
iments  dealt  with  the  effect  of  power  level  changes  as  a  function  of  time,  for  the  seismic  data. 

In  this  experiment  an  exponentially  damped  sinusoid,  whose  frequency  was  1  1 1 z ,  was  generated 

1 1 

at  a  number  of  stations.  The  stations  used  were  the  same  as  those  employed  previously  in 
processing  the  27  October  1%6  Novaya  Zemlya  event.  The  waveforms  at  all  of  the  stations  were 
made  to  be  identical,  indicating  that,  the  data  were  being  generated  by  a  stationary  source  located 
at  a  fixed  point  on  the  earth  which  was  taken  to  be  the  epic  enter  for  the  previously-mentioned 
event.  Thus,  when  the  beamforming  process  is  performed  over  successive  two-second  intervals 
of  data,  the  location  of  the  peak  of  the  output  beam  power  should  remain  stationary  at  the  assumed 
epicenter. 

However,  it  was  found  that  the  location  of  this  peak  did  not  remain  stationary,  but  shifted 
by  a  considerable  amount.  This  shift  of  position  was  due  to  power  level  inequalities  along  the 
waveform.  In  other  words,  it  is  possible  that,  when  two-second  segments  of  the  waveform  at 
each  sensor  are  shifted  using  time  delays  cor  responding  to  a  source  location  other  than  the  ep¬ 
icenter,  the  decrease  in  beam  power  output  due  lo  the  misalignment  of  the  t races  is  more  than 
compensated  for  by  the  increase  in  power  level  in  each  of  the  waveforms.  This  effect,  due  to 
power  level  variations,  is  undesirable  and  should  be  removed,  since  otherwise  there  would  be 
very  little  credence  plac  ed  in  the  source  structure  computation. 

In  order  to  remove  this  effect,  compensation  for  power  level  variations  was  employed  in 
the  beamforming  process.  This  was  accomplished  by  measuring  beam  power  output  as 
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and  I).  ,  is  the  data  in  the  j^1  sensor  at  the  k1*1  time  sample,  t . .  is  the  time  delay  required  in  the 
j  sensor  for  the  i  1  source  location,  N  is  the  number  of  sensors  and  M  is  the  number  of  data 
samples  contained  within  the  integration  time  of  two  seconds,  which  for  20  samples  per  second 
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would  be  equal  to  40.  In  the  previous  beamforming  operation  where  no  power  compensation  was 
used,  the  preceding  formula  applies  with  Wj  .  =  1.  In  essence,  this  new  beamforming  operation 
is  measuring  the  coherence  of  the  delayed  waveforms  and  is  relatively  immune  to  power  level 
variations.  When  this  new  beamforming  operation  was  applied  to  the  artificially  generated  data 
mentioned  previously,  the  location  of  the  peak  power  output  remained  stationary  as  was  desired. 

Another  experiment  was  performed  to  test  the  validity  of  the  P-wave  source  structure  com¬ 
putation  using  the  new  beamforming  operation  which  compensates  for  power  level  variations.  In 
this  experiment  four  source  locations  were  assumed  at  the  points  X,  A,  B,  C,  as  shown  in 
Fig.  II  —  1 .  The  point  X  corresponded  to  the  location  of  the  epicenter  for  the  27  October  1966 
Novaya  Zemlya  event.  A  pulse  is  assumed  to  propagate  from  a  given  source  location  (point  X, 

A,  B,  or  C),  to  each  sensor  in  the  array,  which  is  the  same  as  that  discussed  previously,  with 
a  travel  time  determined  by  the  time  required  for  a  compressional  wave  to  propagate  from  point 
X  to  the  given  source  location  plus  the  time  required  for  a  P-wave  to  travel  from  the  given 
source  location  to  the  sensor.  This  latter  time  is  determined  from  travel  time  tables  by  using 
the  distance  from  the  given  source  location  to  the  sensor.  The  time  required  for  a  compres¬ 
sional  wave  to  propagate  from  point  X  to  a  given  source  location  is  equal  to  the  distance  between 
these  points,  in  km,  divided  by  6  km/sec.  The  amplitudes  of  the  pulses  from  X,  A,  B,  C  were 
assumed  lobe  1,  \r2^/2,  l/2,  \fT/ 4,  respectively,  corresponding  to  successive  3-dB  increments 
in  power  level.  This  sequence  of  four  pulses  at  each  sensor  was  applied  to  a  bandpass  filter, 

0.50  to  1.70  11z,  to  make  the  data  look  somewhat  like  actual  observed  seismic  data. 

The  results  of  the  source  structure  computation  are  shown  in  Fig.  11-1,  which  shows  con¬ 
tours  of  —10  log(P./P  ),  where  P  is  the  maximum  value  of  P..  The  level  of  the  contours 
&  l  max  max  i 

varies  from  0  to  3  dB  in  increments  of  1  dB.  The  dotted  c  ontour  surrounding  the  epicenter, 
marked  X,  in  Fig.  II - 1  indicates  the  locus  of  points  of  possible  sources  of  scattered  P-waves. 
The  results  shown  in  Fig.  II-l  indicate  that  the  location  of  the  peak  power  output  shifts  by  an 
amount  which  is  in  reasonably  close  agreement  with  the  known  shift  in  location  of  the  source  of 
the  P-waves.  That  is,  the  location  of  the  zero  in  Fig.  11-1  shifts  successively  from  point  X  to 
A,  to  B,  and  then  to  C.  Thus,  a  reasonably  accurate  computation  of  source  structure  has  been 
obtained. 

A  more  stringent  test  of  the  method  was  performed  by  assuming  a  shift  of  source  location  in 
two  opposite  directions  simultaneously,  instead  of  a  shift  m  a  single  direction  as  was  done  pre¬ 
viously.  In  this  case,  the  shift  in  the  location  of  the  peak  output  power  was  unable  to  follow  the 
shift  in  location  of  the  source.  Thus,  it  would  appear  that  the  source  structure  computation  is 
useful  only  in  those  situations  where  the  source  location  is  moving  in  a  single  direction.  The 
new  beamforming  operation  has  been  applied  to  several  events,  including  both  earthquakes  and 
presumed  underground  nuclear  explosions.  The  results  will  be  reported  on  soon  in  a  Technical 

Note.  T  _ 

J.  Capon 


B.  SHORT- PERIOD  DATA  EQUALIZATION  AND  DECONVOLUTION 

Measurements  are  being  made  of  large  magnitude,  short-period  teleseismic  P  events  to 
estimate  the  variation  in  crustal  filtering  recorded  at  the  different  subarray  sites  of  LASA.  For 
each  event  crustal  equalization  filters  have  been  computed  to  convert  a  chosen  reference  trace 
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into  various  subarray  traces.  By  repeating  the  calculations  for  other  events  from  the  same  ep- 
icentral  location,  a  similar  set  of  equalization  filters  should  be  obtained.  This  assumes  that 
each  teleseismic  event  has  a  uniform  wave  shape  over  the  wave  front  incident  at  the  base  of  the 
crust  and  that  the  recorded  signal  variations  over  the  aperture  of  1  ASA  are  caused  predominantly 
by  variations  in  crustal  structure  under  LASA.  This  appears  to  be  true  from  the  data  studied 
which  is  discussed  below. 

Eight  Aleutian  events  of  magnitude  5.4  or  greater  were  selected  for*  the  initial  study  of 
crustal  filtering.  These  events  are  listed  in  Table  11  —  4  in  order  of  increasing  epicentral  distance 
from  LASA.  Events  1  and  2  are  from  the  Fox  islands;  Events  3  through  6  are  closely  located 


TABLE  11-1 

ALEUTIAN  ISLANDS  P  EVENTS  RECORDED 

AT  LASA 

Event 

Date 

Arrival  Time 

Distance 

Back  Azimuth 

Phase  Velocity 

Magnitude 

1 

12/4/65 

2:19:39  GMT 

39.6° 

299.7° 

13. 39  km/sec 

5.7 

2 

1/14/68 

12:48:34 

42.9° 

303. 1° 

13.66 

5.6 

3 

6/22/67 

15:44:49 

45.3° 

304.1° 

13.94 

6.3 

4 

2/25/68 

18:16:30 

45.5° 

302.4° 

13.96 

5.6 

5 

12/19/67 

14:48:52 

46.0° 

303.9° 

14.03 

5.5 

6 

9/16/67 

8:40:06 

46.9° 

303.7° 

14. 14 

5.8 

7 

9/13/67 

18:50:10 

51.2° 

308.9° 

14.71 

5.6 

8 

2/9/68 

15:42:10 

52.2° 

311.5° 

14.85 

5.4 

in  the  Andreanof  Islands,  and  Events  7  and  8  are  from  the  Near  Islands.  The  clearest  variations 
of  crustal  filtering  in  this  set  of  events  are  seen  on  the  subarray  traces  for  Event  3  which  are 
shown  in  Fig.  11-2.  The  variations  of  amplitude  and  time  duration  of  the  P-wave  shown  here  are 
similar  to  variations  observed  for  the  other  events  in  this  group. 

In  order  to  calculate  equalization  filters  for  the  subarrays  of  LASA,  each  recorded  event  is 
assumed  to  have  the  form 

e^t)  =  c.(t)*  s(t) 

where  s ( t )  is  the  incident  teleseismic  plane  wave  source  at  the  base  of  the  crust  under  LASA, 
c.(t)  is  the  crustal  transfer  function  for  subarray  site  i,  and  *  denotes  convolution.  The  crustal 
transfer  function  may  depend  upon  azimuth  and  distance  of  the  event.  An  equalization  filter  for 
converting  e^(t)  to  e^t)  for  an  event  is  given  in  the  frequency  domain  by 

Fd2(w)  R2(a.')/Ei(o i)/C^u) 

where  a;  is  frequency  and  (E.(a>),  C.(o>)]  are  Fourier  transforms  of  [e.(t),  c . (t )].  Note  that  the 
equalization  filter  is  independent  of  the  source  waveform  if  the  souree  location  is  fixed. 

Equalization  filters  were  computed  for  all  subarray  sites  for  each  of  the  eight  events.  The 
reference  trace  used  was  the  sum  of  the  134  and  AO  traces.  This  trace  generally  had  a  simple 
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form  and  was  used  for  this  reason  as  a  reference  trace.  The  reference  traces  for  the  eight 
events  are  shown  in  Fig.  II - 3 ,  plotted  against  the  epicentral  distance  of  each  event  from  LASA. 

The  equalization  filters  were  computed  in  the  time  domain  by  a  least-squares  technique. 

Two  examples  of  the  filter  computation  are  shown  in  Figs,  ll-4(a)  and  (b).  In  (a)  a  filter  is  cal¬ 
culated  to  shape  five  seconds  of  the  input  trace  of  Event  3  into  10  seconds  of  a  ringing  P-wave 
recorded  at  site  Bl.  This  equalization  filter  is  applied  to  the  complete  input  trace  and  the  re¬ 
sult  compared  with  the  complete  desired  output  trace.  It  is  of  interest  that  the  equalization 
filter  equalizes  the  pP  phases  as  well  as  the  P  phases  on  the  two  traces.  This  is  another  con¬ 
sistency  check  on  the  data  since  the  pP  and  P  phases  traverse  almost  identical  paths  to  the 
seismometers  in  LASA.  In  Fig.  II-4(b)  a  more  complex  source  is  apparent  on  the  input  trace  of 
Event  4.  The  equalization  procedure  again  equalizes  later  portions  of  the  traces  not  used  to  de¬ 
sign  the  filter. 

Finally,  it  is  useful  to  examine  the  similarities  of  equalization  filters  for  different  events. 
Figure  11-5  shows  equalization  filters  designed  to  shape  the  AO  plus  B4  reference  trace  to  three 
ringing  subarray  sites,  i.e.,  (a)  Bl,  (b)  C2,  and  (c)  Dl.  The  vertical  scale  for  the  filters  equals 
1  between  the  pips.  There  are  several  interesting  features  of  these  filters  to  notice.  A  strong 
multiple  spike  appears  on  all  three  filters  for  the  closely-spaced  events  at  A  ~  46°.  The  equal¬ 
ization  filters  for  C2  have  at  least  two  reliable  multiple  spikes  for  these  events  even  though  the 
input  traces  shown  in  Fig.  II - 3  are  quite  different  in  shape.  Events  7  and  8,  which  are  only  5° 

away  from  Events  3  through  6  seem  to  lose  the  multiple  spikes.  This  verifies  the  extreme  sen- 

12 

sitivity  of  the  crustal  response  to  azimuth  and  phase  velocity  changes  as  noted  by  Mack  and 

1  3 

Earner  for  nonuniform  crustal  structures. 

Results  of  this  preliminary  work  show  that  crustal  equalization  filters  can  be  designed  for 
teleseismic  events  from  localized  regions.  A  possibility  now  being  investigated  is  the  deconvo¬ 
lution  of  all  subarray  traces  back  to  a  short  duration  signal.  This  inverse  equalization  filtering 
may  be  useful  for  short -period  detection  of  pP  phases  buried  in  P-wave  codas.  In  the  light  of 
the  strong  crustal  filtering  variations  shown  above,  it  is  clear  that  short-period  data  can  be  used 
for  discrimination  only  with  the  understanding  that  the  discriminant  may  be  site-dependent. 

C.  W.  Frasier 
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LONGITUDE  (Host) 


(c)  A  to  6  SEC 


(d)  6  to  8  SEC 


Fig.  11-1. 


P-wave  source  structure  result  far  artifically-generated  seismic  data. 
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TIME  (s«c) 


Fig.  11-2.  Summed  subarray  traces  at  LASA  far  Event  3,  a  magnitude  6.3  P-wave 
from  the  Andreanof  Islands. 
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Fig.  11-3.  Aleutian  Islands  events  used  ta  study 
the  crustal  filtering  af  P-waves  at  LASA.  Each 
trace  is  the  sum  af  subarray  sites  B4  and  AO. 
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Fig.  11-4.  Examples  af  equalization  procedure  for 
LASA  subarray  sites  far  twa  Aleutian  P  events;  in 
each  case  5  sec  af  the  (B4  +■  AO)  trace  sum  is  used 
as  the  input  far  the  filter  calculation,  (a)  Event  3, 
equalization  ta  B1  subarray,  (b)  Event  4,  equali¬ 
zation  ta  El  subarray. 
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-2-9063 


-  \z"\AA. - 


(c)  D1 


Fig.  11-5.  Examples  af  equalization  filter  impulse  responses.  In  each  case 
the  input  is  the  B4  +  AO  subarray  sums  af  Fig.  11-2.  Desired  output  subarrays 
are  (a)  Bl,  (b)  C2,  and  (c)  Dl. 
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III.  MISCELLANEOUS  PROJECTS 


A.  VESPA  AS  A  GEOPHYSICAL  TOOL 

A  study  has  been  made  to  determine  if  the  program  VESPA  (Velocity  Spectral  Analysis)  will 

be  of  use  in  measuring  rav  parameters  at  LASA  of  phases  whose  relative  arrival  times  cannot  be 

measured  accurately  by  eye  because  of  single  trace  amplitudes  that  are  too  weak  or  near  other 

14  l  5 

arrivals.  VESPA  is  described  in  the  June  1968  and  December  1968  SATS.  ’  It  steers  the 
array  to  the  appropriate  azimuth  and  forms  61  beams  with  slowness  from  0.0  to  14.0  sec/deg 
in  increments  of  0.4  sec/deg.  The  powers  of  each  of  the  beams  are  computed  over  1-sec  inter¬ 
vals  and  these  powers  are  contoured  as  a  function  of  time  and  slowness.  Using  the  subarray 
outputs  of  the  D,  E,  and  F  rings  of  LASA,  the  scheme  has  been  applied  to  the  P  arrivals  from 
17  earthquakes  within  the  distance  range  of  140°  to  170  .  1'rom  events  within  this  range  multiple 
arrivals  are  observed  due  to  the  earth's  core  and  its  internal  structure. 

As  an  example.  Fig.  1 1 1  —  1  shows  a  contour  plot  of  power  vs  slowness  and  time  (hereafter 
calked  a  vespagram)  for  an  earthquake  of  October  45,  1967.  The  USGGS  E DR  parameters  for 
this  event  are:  origin  time  i5h  09m  3  5.3  s,  location  50.15r,S  and  114.47°E,  and  normal  depth 
(restricted  to  33  km).  With  respect  to  the  LASA  centerpoint  the  event  is  at  a  distance  of  153.1 
and  an  azimuth  of  44  7.7  .  The  beam  formed  using  the  D,  E,  and  F  subarray  sums  delayed  ar¬ 
bitrarily  for  a  slowness  of  4.8  sec/deg  is  shown  beneath  the  vespagram  in  Fig.  111-1.  T  hree 
distinct  arrivals  are  apparent  and  numbered  on  the  beam  trace.  The  vespagram  computed  from 
this  data  and  contoured  in  4-dB  increments  down  from  the  peak  power  shows  three  peaks  asso¬ 
ciated  with  these  arrivals,  each  at  a  different  slowness.  The  time  intervals  between  peaks  1 
and  4,  and  1  and  3  are  about  7  and  40  sec,  respectively.  The  core  phase  travel  time  tables  of 
Bolt1  give  arrival  times  for  three  phases  at  153°.  The  earliest  and  latest  of  these  are  bast'd 
on  observations,  the  intermediate  phase  is  a  smooth  extrapolation  of  observations  at  lesser 
distances  and  is  based  on  a  theoretical  model.  At  153°  the  time  intervals  between  the  first  and 
second  and  first  and  third  phases  of  Bolt's  tables  are  4.8  and  18.1  sec.  At  154  these  intervals 
are  6.0  and  41.0  sec.  Slowness  values  computed  from  Bolt's  tables  at  153°  are  1.4  5,  4.4,  and 
4.45  sec/deg.  For  the  event  shown  in  Fig.  Ill  —  1 ,  values  of  1.4,  4.6,  and  4.45  sec/deg  were 
measured  on  a  vespagram.  In  Fig.  1 1 1  —  1  the  contouring  interval  is  greater  than  that  normally 
used  in  order  to  make  the  display  clearer  in  reproduction. 

The  vespagrams  of  16  additional  events  were  "read"  in  a  similar  manner.  The  measured 
slownesses  are  plotted  (open  circles)  against  measured  travel  time  in  Fig.  II 1-4  in  order  to  con¬ 
voy  as  much  information  as  possible  concerning  the  relevancy  of  the  measurements.  FSCGS 
data  was  used  for  origin  times  and  locations;  all  of  the  depths  were  either  restricted  to  or  com¬ 
puted  to  be  less  than  3  3  km.  Peaks  in  the  velocity  spectrum  were  timed.  These  usually  ranged 
from  1  to  4  sec  following  the  onset  of  the  group  of  contours  indicating  the  arrivals.  Curves 
computed  from  Bolt's  tables  arc  drawn  on  1  ig.  111-2  as  a  standard  against  which  to  compare  the 
data.  The  letters  refer  to  branches  given  in  his  tables.  The  curves  are  dashed  where  his  tables 
indicate  extrapolations. 

it  must  be  stressed  that  the  plotted  points  of  Fig.  III-4  represent  a  rough  test  of  an  analysis 
scheme  rather  than  a  refined  set  of  observations.  The  plotted  observations  follow  the  established 
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curve  of  PKP^  well  at  larger  slownesses.  At  smaller  slownesses,  along  the  PK1KP  branch, 
the  analysis  appears  to  trend  toward  a  greater  value  for  d/dT(dT/dA)  than  indicated  by  the 
established  observations. 

In  conclusion,  the  technique  appears  promising  as  a  geophysical  tool.  However,  the  assign¬ 
ment  of  error  ranges  and,  in  the  present  case,  station  corrections  at  the  higher  velocities  re¬ 
main  as  problems.  .  . 


B.  ANALYSIS  OF  MICROBAROGRAPH  DATA  FROM  ATMOSPHERIC  EXPLOSIONS 

The  French  thermonuclear  device,  exploded  on  24  August  1968  in  the  South  Pacific,  was 

recorded  by  the  Large  Aperture  Mierobarograph  Array  (LAMA)  situated  in  Montana.  High- 

9 

resolution  frequency-wavenumber  techniques  have  been  applied  to  the  recordings  and  some 

examples  are  given  in  Figs.  III-3(a)  and  III  —  3(b).  The  two  frequencies  represented  correspond 

1 7 

to  maxima  in  the  frequency  spectrum  of  the  received  signal.  The  phase  velocity  of  the  0.006 -Hz 
wave  is  seen  to  be  slightly  over  0.3  km/sec.  This  agrees  with  the  accepted  value  for  acoustic 
propagation  in  the  lower  atmosphere.  The  velocity  for  the  0.01-IIz  wave  is  apparently  somewhat 

hlghcr'  H.  Mack  (M.I.T.  Earth 

and  Planetary  Sciences  Department) 


C.  PDP-7  DATA  PROCESSING  FACILITIES  DEVELOPMENT 

The  following  hardware  items  have  been  added  to  our  computer  system: 

(1)  The  Beta  hard  copy  printer  has  been  delivered  and  installed  so  that  it 
can  be  shared  by  the  two  computers.  In  a  few  seconds  this  device  can 
generate  8-1/2-  x  11 -inch  copies  of  any  data  which  has  been  displayed 
on  our  computer-driven  D.E.C.  340  scopes. 

(2)  The  40.8-kilobaud  data  link  between  the  M.I.T.  Information  Processing 
Center  (IPC)  and  one  of  the  PDP-7's  has  been  installed  and  preliminary 
closed  loop  tests  that  exercise  the  PDP-7  interface,  phone  line  and  the 
tw'o  modems  show'  satisfactory  performance.  The  IBM  360/40  at  the 
IPC  is  being  modified  to  accept  data  in  our  format  and  should  soon  be 
ready  for  remote  job  entry  of  large  data  files. 

(3)  A  shaft  encoder  and  knob  box  have  been  constructed  and  installed  on 
computer  2,  along  w'ith  a  32-ehanncl  multiplexer  and  A/D  converter. 

As  a  result,  both  computers  can  use  the  Data  Analysis  Console  at  the 
same  time  to  analyze  different  data. 

(4)  The  design  modifications  for  the  3:1  speedup  for  the  D.E.C.  340  display 
have  been  completed.  It  was  found  to  reduce  the  flicker  by  at  least  a 
factor  of  three,  and  can  be  increased  to  five  at  the  expense  of  decreasing 
the  point  plotting  accuracy  slightly. 

The  following  items  of  software  have  been  completed  or  are  in  progress: 

(5)  The  New  Data  Analysis  Console  has  been  designed  and  several  sub¬ 
programs  have  been  coded.  The  drum  has  been  returned  to  the  factory 
for  replacement  of  the  spinner  and  flying  head  assembly  to  correct  a 
manufacturing  error.  Consequently,  final  checkout  and  debugging  is 
deferred  until  the  drum  is  back  in  operation. 

(6)  The  software  for  the  data  link  to  the  360  is  complete  so  that,  when  the 
IBM  360  hardware  is  ready,  we  will  be  able  to  send  jobs  to  the  360  on 
a  time-shared  basis  with  the  New  Data  Analysis  Console. 
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(7)  A  program  has  been  written  that  makes  hard  copy  prints,  via  the  Beta 
device  and  PDP-7  computer,  from  magnetic  tapes  made  by  the  IBM  360 
in  Stromberg-Carlson  4020  format  for  plotting  by  the  latter  device. 

This  program  has  been  included  in  the  package  called  "Routine  Processing" 
as  described  in  the  previous  SATS.*8 

(8)  The  FORTRAN  compiler  tape  backup  system  has  been  completed  and  is 
now  usable.  The  final  system  will  be  put  on  the  drum  and  should  be  easier 
to  use  and  faster.  The  IBM  G  level  FORTRAN  syntax  checker  is  just 
being  started  by  the  subcontractor.  The  compiler  has  been  documented 

by  us  and  is  published  as  Chapter  IV  in  the  PDP-7  Software  Handbook. 
Several  FORTRAN  user  packages  are  currently  being  programmed  in 
machine  language  to  give  our  system  capabilities  similar  to  software 
routines  in  the  IBM  360  system.  Those  include  Stromberg-Carlson-type 
plotting  routines  and  Input/Output  routines  to  read  and  write  LASA  format 
data  tapes.  This  will  enable  many  of  the  FORTRAN  programs  which  are 
currently  being  run  on  the  IBM  360  to  be  run  on  our  system.  This  should 
encourage  the  scientist -type  programmers  to  use  our  system  instead  of  the 
IBM  360  because  they  will  be  able  to  read  data  tapes  and  create'  displays 
without  having  any  machine  language  experience.  This  has  the  added  ad¬ 
vantage  of  instantaneous  output  which  should  result  in  more  efficient 

utilization  of  human  resources.  ..  ...  , 

P.  Meek 


D.  SEISMOLOGICAL  SUBROUTINE  LIBRARY 

A  library  has  been  established  that  consists  entirely  of  subroutines  useful  for  general  seis- 
mologieal  calculations.  Because  the  setup  and  structure  of  this  library  is  very  similar  to  a 
computer  operating  system  library  it  will  allow  the  programmers  rapid  and  easy  access  to  a 
variety  of  subroutines  required  in  seismologieal  problems.  The  structure  of  most  of  the  sub¬ 
routines  is  rather  flexible  and  at  the  same  time  permits  them  to  be  relatively  free  of  input  and 
output  operations.  Because  of  the  variety  of  subroutines  that  are  available,  a  great  amount  of 
programmer  time  can  be  saved  by  the  use  of  these  self-contained  and  tested  subroutines. 

Included  in  the  variety  of  subroutines  that  are  available  in  the  seismologieal  subroutine 
library  are  the  basic  routines  to  compute  distance,  azimuth,  travel  time  of  body  waves,  magni¬ 
tude,  depth  of  the  hypocenter  based  on  pP,  elevation  station  correction,  ellipticity  corrections, 
Montana  LASA  time  station  corrections  and  the  name  of  the  geographic  region  at  any  given  lat¬ 
itude  and  longitude.  Other  subroutines  are  useful  with  data  from  large  arrays.  These  routines 
will  allow  the  programmer  to  obtain  azimuth,  horizontal  phase  velocity  and  distance  to  the  epi¬ 
center  based  on  a  least-squares  fit  of  a  plane  wave  to  the  observed  arrival  times.  The  plane 
wave  delays  across  an  array  can  also  be  obtained  from  any  velocity  and  azimuth.  Another  basic 
routine  will  eonvert  an  array  of  station  positions  of  latitude  and  longitude  into  positions  in  an  xy 
plane.  A  general  fast  Fourier  transform  routine  is  also  available  for  the  purpose  of  obtaining 
s  pe  ct  ra . 

There  are  several  routines  available  to  aid  the  programmer  with  plotting  and  interpretation 
of  data.  One  of  the  most  useful  routines  will  develop  contour  maps  of  data  sampled  in  xy  coor¬ 
dinates  or  in  polar  coordinates.  The  program  generates  a  tape  which  can  he  plotted  on  the 
Stromberg-Carlson  4020  plotter  or  by  the  Beta  attached  to  the  PDP-7  system.  As  a  further 
aid  to  interpreting  the  contour  plot,  there  is  a  routine  that  will  plot  the  two-dimensional  array 
as  a  three-dimensional  figure.  Other  subroutines  will  allow  the  programmer  to  plot  waveforms, 
long-period  or  short -period,  and  specify  any  time  scale  desired.  Several  other  routines  are 
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available  that  perform  basic  functions,  such  as  reading  or  writing  data  tapes  recorded  at  several 
large  arrays  or  plotting  data  on  the  printer  for  a  quick  look  at  functions  or  data. 

Because  of  the  variety  of  programs  in  the  library  and  their  use  in  basic  seismologieal  prob¬ 
lems,  it  is  hoped  that  the  programmer  can  be  relieved  of  the  effort  required  to  write  his  own 
basic  routines.  The  library  is  permanently  mounted  in  the  M.I.T.  Information  Processing 
Center  operating  system  and  is  usable  by  programmers  utilizing  the  IBM  360  at  the  IPC. 

R.  M.  Sheppard 
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1 1 1-1.  Beam  farmed  using  the  D,  E,  and  F  ring  subarray  sums  delayed  at  2.8  sec/deg  and  the  vespagram 
puted  from  the  same  data  far  the  event  af  25  October  1967. 
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TRAVEL  TIME 


Fig.  II 1-2.  Slowness  versus  travel  times  read  from  vespagrams  of  core  phoses 
ot  LASA  from  1  7  events  at  140°  <  A  <  170°. 
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Fig.  111-3.  High  resolution  wavenumber  spectrum  of  acoustic  waves  generated  by  French  device 
on  24  August  1968  and  recorded  at  the  LAMA,  (a)  0.006  Hz.  (b)  0.01  Hz. 
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